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ABSTRACT 

The  use  of  low  thrust  devices  with  continuous  and  discontinu- 
ous thrust  programs  is  investigated  to  determine  whether  or  not  a 
discontinuous  thrust  program  will  provide  a  greater  payload  and 
structure  mass  fraction,    where  the  program  must  impart  a  specified 
amount  of  energy  to  the  satellite  in  a  specified  total  time,    starting 
from  a  circular  parking  orbit.     A  discontinuous  thrust  program  is 
developed,   based  on  the  elliptic  orbit,    using  a  perturbation  analysis, 
and  a  series  solution  is  obtained  which  permits  investigation  of  this 
program  up  to  angles  of  one  radian  either  side  of  the  perigee  of  the 
elliptic  orbit.      These  results  are  compared  with  a  continuous  thrust 
program  which  gave  a  spiral  orbit.     Under  certain  conditions,    where 
storage  batteries  must  be  carried  as  part  of  the  payload  but  are 
available  for  use  during  the  thrust  program,    the  discontinuous  thrust 
program  is  found  to  provide  a  greater  payload  and  structure  mass 
fraction.     Further  investigation  at  angles  greater  than  one  radian 
either  side  of  perigee  appears  to  be  warranted. 
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g  gravity  at  surface  of  central  mass 

R  radius  at  surface  of  central  mass 
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r  radius  at  any  point  in  the  orbit 

0  angle  of  position  vector  in  the  plane  of  motion 

P  angle  of  velocity  vector  in  the  plane  of  motion 
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h  non-dimensional  angular  momentum 
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m,  M  mass 

M  initial  mass  of  satellite 
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M  mass  of  power  supply 
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P 
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IvL  mass  of  payload 
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h/L  mass  of  batteries 

K,  specific  mass  of  batteries 

K  specific  mass  of  power  supply 
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(    )  vector  notation 
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(    )  derivative  with  respect  to  time 

PR  perigee  restriction 

\  Lagrange  multiplier 

r\  specific  power 
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N  number  of  thrust  cycles 

2 

FR   _  factored  ratio  of  v 

Y 

m  ratio  of  M,  to  M 

^  oatt.  g 

a  energy  storage  factor 


I.     INTRODUCTION  AND  SUMMARY 

Once  in  a  parking  orbit,   low  thrust  devices  may  be  used  to 
alter  the  path  of  a  satellite     In  particular,    it  is  anticipated  that  one 
use  of  the  low  thrust  device  will  be  to  add  energy  to  the  satellite  for 
an  escape  or  a  near  escape  condition.      The  exact  path  and  thrust 
program  must  be  determined  by  the  mission  requirement,    but  in  any 
event,   the  system  should  be  optimized  for  efficient  launching. 

The  problem  of  optimizing  an  arbitrary  mission  is  a  for- 
midable task  to  say  the  least.      To  simplify  matters,    it  is  assumed 
that  the  mission  is  to  impart  a  specified  energy  in  a  specified  total 
time  to  the  maximum  fraction  of  payload  and  structure  mass.      To  do 
this,    one  must  make  efficient  use  of  the  power  supply  and  propellant. 

The  usual  variation  procedures  may  be  used  to  derive  the 
differential  equations  for  optimum  direction  and  magnitude  of  thrust 
when  continuous  thrust  is  employed,    as  shown  by  Irving  in  reference 
(1).     He  has  shown  that  a  constant  acceleration  is  optimum  for  the 
gravity-free  case.     No  analytic  solution  has  been  found  for  the  dif- 
ferential equations  in  the  central  force  field,   but  Irving  has  made  a 
numerical  study  and  found  the  acceleration  to  be  practically  constant 
for  various  escape  programs. 

In  reference  (2),    Casey  has  shown  the  distinct  advantage  of 
energy  addition  by  thrusting  near  the  perigee  of  an  elliptic  orbit  in 
such  a  manner  that  the  perigee  remains  constant.     His  thrust  pro- 
gram is  formulated  to  require  that  the  perigee  distance  remain  con- 
stant throughout  the  thrust  cycle,    and  results  are  obtained  by 
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numerical  procedures.      The  reason  for  thrusting  near  the  perigee  is 
to  take  advantage  of  the  greater  velocity,    since  in  any  elliptic  orbit 
the  velocity  is  greatest  at  the  perigee  and  since  the  energy  addition 
rate  is  proportional  to  the  velocity.     Further,    if  the  perigee  distance 
is  kept  constant,   then  the  perigee  velocity  will  increase  as  the  eccen- 
tricity of  the  elliptic  orbit  increases. 

As  pointed  out  by  Casey,    there  will  be  a  period  of  coasting 
between  each  thrust  cycle  which  will  make  a  considerable  contribu- 
tion to  the  total  time.     For  this  reason,    we  cannot  say  simply  by  in- 
spection whether  the  continuous  or  discontinuous  method  of  energy 
addition  will  provide  a  greater  fraction  of  payload  and  structure 
mass,    where  both  methods  are  required  to  impart  a  specified  energy 
in  a  specified  total  time. 

To  investigate  this  problem,   we  shall  use  the  idea  of  re- 
stricting the  perigee  distance,   but  only  at  the  beginning  and  end  of 
each  thrust  cycle.      This  permits  added  freedom  in  an  effort  to  gain 
more  energy.     Since  we  are  only  considering  low  thrust  devices,    a 
perturbation  technique  is  used  with  a  simple  variation  to  develop  the 
first  perturbation  integral  equations  for  the  maximum  change  of  en- 
ergy during  one  thrust  cycle  and  for  the  Lagrange  multiplier  neces- 
sary to  satisfy  the  above  perigee  restriction.      The  equations  apply 
for  any  period  of  thrusting  up  to  one  complete  circuit  of   2tt   radians, 
but  could  not  be  integrated  in  closed  form.     By  making  a  series  ex- 
pansion of  the  integrands,    the  equations  may  be  integrated  term  by 
term  to  obtain  results  useful  to  an  angle  of  one  radian  either  side  of 
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the  perigee.      While  this  does  not  permit  a  complete  comparison  of 
these  two  methods  of  energy  addition,    it  is  possible  to  draw  definite 
conclusions  for  angles  less  than  one  radian. 

After  obtaining  equations  for  one  cycle  of  energy  addition, 
the  effect  of  a  large  number  of  cycles  is  obtained  by  integrating  the 
changes  of  the  various  elliptic  orbit  parameters. 

To  compare  these  results,    an  approximate  solution  for  the 
spiral  orbit,   using  a  low  thrust  device,   is  obtained.     Then  a  compari- 
son is  made  of  the  payload  and  structure  mass  fraction,    where  both 
methods  are  required  to  fulfill  the  mission  of  imparting  a  specified 
energy  in  a  specified  total  time. 

Finally,    since  the  discontinuous  thrust  program  could 
generate  energy  while  coasting,  the  use  of  storage  batteries  to  aug- 
ment the  power  available  during  the  thrust  cycle  will  be  considered. 
In  one  case,    the  batteries  will  be  considered  as  part  of  the  power 
supply;  and  in  another,    they  will  be  taken  as  part  of  the  payload,   but 
available  for  use  during  the  thrusting  program. 

When  thrust  is  only  permitted  up  to  an  angle  of  one  radian 
either  side  of  perigee,    the  discontinuous  thrust  program  was  found 
to  have  no  payload  and  structure  mass  fraction  advantage  unless  bat- 
teries are  employed.      When  the  silver-zinc  battery  is  considered  as 
part  of  the  power  supply,    there  is  an  area  of  marginal  advantage  due 
to  the  low  specific  mass  of  this  battery,   but  its  poor  recharging  re- 
liability rules  out  its  present  use  in  such  a  system. 

The  nickel-cadmium  storage  battery,    which  is  very  relia- 
ble,  provided  an  area  of  definite  advantage  for  the  discontinuous 
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thrust  program  when  the  battery  is  considered  as  part  of  the  payload. 
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II.     PRELIMINARY  ANALYSIS 

We  shall  first  write  the  equations  of  motion  for  the  satellite, 
making  the  usual  simplifying  assumptions  of  a  very  large,  central, 
spherical  mass  and  take  the  coordinate  system  at  the  center.     We 
shall  further  assume  motion  in  a  plane,    keeping  the  thrust  in  the 
plane  of  the  orbit,    and  neglect  the  effects  of  other  bodies.     Then, 
using  plane  polar  coordinates  ,    as  in  fig.    1,   and  the  dot  notation  for 
time  derivative,  we  have: 

n„   =  7  n  (2.i) 

/f>*    =  Th     +■  j  ne  <2-2) 


n, 


.  T(ri-ndz)  +  j^^Jn'e])  (2.3) 


Equating  equation  2.  3  to  the  acceleration  due  to  thrust,  q,  »    and  grav- 
ity,  we  have: 

+    j(a.    5<N(f34-d>))  (2.4) 

where     RQ    is  the  radius  of  the  central  mass,    and    <%0     is  the  acceler- 
ation of  gravity  at    n  =  R0   .     Since     a/%°    is  small  for  a  low  thrust 
device,    we  shall  first  take     0.  -  o         and  obtain  the  usual  solution  in 
dimensionless  form.      The  thrust  will  be  considered  later  as  a  per- 
turbation.    Choosing     &     as  the  independent  variable,    and  letting  (  )' 
be  differentiation  with  respect  to    &    ,    then  with    U     as  the  dependent 
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variable  and    h     as  the  dimensionless  angular  momentum,    we  h 


ave: 


u     =     7T  (2-5) 


nze  ~-  )\URl  f\o  ft 


(2.6) 


but  when     a  =  O      *    jz.(  J)l&  )  =  0  ',  therefore     h  =  constant  = 


"Mo    • 


©   =  l/i7i?7  f»e  u  (2.7) 


ri      =  -]/?«Rj-Pl,  l/  (2.8) 


h'    =  -  3.  ft.  tr  u 


/76^=     3ot\*U: 


(2.9) 


(2.  10) 


Substituting  into  equation  2.4,    taking      (X  -  O         ,    and  simplifying, 
we  have: 

U"    +■     U      =     "£2  (2.11) 

A  solution  may  be  written 

U   =     U*  =  ~tfll   +  6°  cos(0-0p)j  (2.12) 

where       c?0      is  the  eccentricity  and       &{>     is  the  angle  of  the  perigee 
position.      The  non-dimensional  time,     A~V     ,   for  any  portion  of  the 
orbit  is  given  by: 
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fa 


(2.13) 


using  equations  2.  7  and  2.  12.      The  integral  can  be  evaluated  to  give: 


In  equations  2.  14  through  2.  26,  we  shall  essentially  follow  the 
work  of  Irving,  reference  (1),  since  his  approach  provides  a  conveni- 
ent parameter  for  computing  the  mass  fraction  of  the  satellite. 

The  mass  of  the  power  supply,    Mw     ,    is  assumed  to  be  pro- 
portional to  the  power  supplied,     P      : 

Mw  =  KP  <2-14) 

where      \(      is  the  specific  mass  of  the  power  supply.     It  is  further 
assumed  that  the  maximum  power  is  utilized  when  thrusting,    appear- 
ing as  kinetic  energy  in  the  exhaust: 

P     =^(-v^)C2  (2.15) 

where  C  is  the  exhaust  velocity,  vv  the  mass  of  the  vehicle  at 
time  +  ,  and  -vvv  is  the  propellant  flow  rate.  We  also  have  the 
acceleration  of  the  satellite  due  to  reaction: 

thrust  YKI  C 


a    = 

mass  yvi 

Using  equations  2.  15  and  2.  16,    we  eliminate       C     : 


(2.  16) 
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Equation  2.  17  may  be  integrated  to  obtain  an  expression  for  mass  as 
a  function  of  time: 

i  I       4-    >    C  G'    At 

TT=7T0  +  I  )  T  (2-18) 

o 

where  M0  is  the  initial  mass  of  the  satellite.  If  we  let  Mp  be 
the  propellant  mass,  then  the  "burn  out"  mass,  Mly  ,  at  time  7£ 
is  given  by: 

M4=M.-M,  (2-19) 

Substituting  in  equation  2.  18  and  using  2„  14,   we  have: 

Mfc  2     Mw 

To  optimize  the  power  supply,   we  write  the  relation  between 
the  payload  mass,      M,       ,    the  structure  mass,      M$       ,    and  the  other 
masses: 

M,  +MS+  Mp+M„  =  M0  (2.21) 

Further,  let 

Tb 
y*  =     Ji  \  aa  Ji*  (2.22) 

o 
Then,    manipulating  equations  2.19,    2.20,    2.21,   and  2.22,   we  have: 


Mo  N.    [(Mw/Mj  +  Y1 


(2.23) 
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There  is  an  infinite  number  of  programs  for     ~ci    which  will  keep    T2 
constant  while  fulfilling  the  mission  as  Mw/M0        is  varied.      Thus, 

taking  the  derivative  of  equation  2.  23  equal  to  zero,   we  have  the 
maximum  of    C  M ,  +  M*)/  M «  with  respect  to       IM^/Mfl      when: 

^=Y-Y^  (2.24) 

Substituting  equation  2.  24  into  2.  23,    we  have: 

-'  -  M$-  =    (l-r)2  (2.25) 

Mo 

Using  equations  2.  21,    2.  24,    and  2.  25,   we  have: 

Hl=  r  (2. 26) 

From  equation  2.  26  we  see  the  range  of      Y     must  be  zero  to  one, 
since      Mp/Mo  must  lie  in  this  range.     Further,   for  a  particular 

mission,   we  see  from  equation  2.  25  that  the  structure  and  payload 
mass  fraction  will  be  greatest  when     Y      is  a  minimum.     After  de- 
termining the  minimum  value  for     Y     ,   the  optimum  power  supply 
mass  is  given  by  equation  2.  24. 

We  are  now  ready  to  investigate       equation  2.  22.     In  reference 
(1),   Irving  has  shown  analytically  that  the  acceleration  must  vary 
linearly  with  time  in  gravity-free  space,   if     V     is  to  be  a  minimum. 
Further,   when  the  mission  of  specified  energy  in  a  specified  time  is 
considered,    the  optimum  acceleration  becomes  constant.     For  the 
case  of  a  central  force  field,    Irving  made  a  numerical  study  using 
continuous  thrust,    and  found  the  optimum  acceleration  to  be  practi- 
cally constant.      Thus,    the  spiral  orbit  with  constant  acceleration  is  a 
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"near  optimum"  solution  which  will  provide  a  fair  comparison  be- 
tween continuous  and  discontinuous  thrust  programming. 

If  we  consider  a  discontinuous  thrust  program  applied  to  an 
elliptic  orbit,   it  seems  clear  that  the  thrust  should  be  applied  near 
the  perigee,    where  the  velocity  is  greatest,    to  impart  the  greatest 
energy,,     If  the  thrust  is  parallel  to  the  velocity  vector,   we  will  im- 
part the  maximum  energy  during  any  given  thrust  cycle.     Such  a  pro- 
gram will  allow  the  perigee  distance  to  increase,  and  with  this  in- 
crease there  is  a  decrease  of  perigee  velocity,   thus  decreasing  the 
energy  imparted  during  a  subsequent  given  thrust  period.      To  obtain 
the  advantage  of  thrusting  near  the  perigee,    but  still  keeping  the 
same  distance  for  a  velocity  advantage  in  subsequent  cycles,    a  pro- 
gram must  be  found  that  imparts  maximum  energy  during  any  given 
thrust  cycle  while  providing  no  net  change  in  the  perigee  distance. 
After  obtaining  the  restricted  perigee  equations  for  one  cycle,   the 
programming  of  thrust  time  will  be  considered  and       Y    will  be  com- 
puted for  a  large  number  of  cycles. 

After  developing  the  restricted  perigee  equations,   we  shall 
obtain  an  approximate  solution  for  the  spiral  orbit  with  continuous 
thrust  to  compare  with  the  restricted  perigee  thrust  program.     Fi- 
nally,  we  shall  consider  the  use  of  batteries  with  the  restricted  peri- 
gee thrust  program.      The  batteries  can  be  charged  while  coasting  to 
provide  added  power  when  thrusting.     In  one  case,    the  batteries  will 
be  considered  as  part  of  the  power  supply,    and  in  another  case  they 
will  be  considered  as  part  of  the  payload. 
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III.     PERTURBATION  EQUATIONS 

When     (X     is  not  zero,   we  have  a  small  non-dimensional  ac- 
celeration: 

£=  ~   <<  I  (3-D 

since  we  are  considering  only  low  thrust  devices.     This  suggests  the 

use  of  a  perturbation  analysis  for  each  cycle. 

Since     7?      is  no  longer  constant,   we  must  rewrite  equation 

2.9: 

H    =    -ji^rtVt^")  (3.Z) 

and  for  the    J       component  of  equation  2.  3  we  also  have: 

JL  ±lnze)  =  g.u'AA'  (3-3> 

By  using  fig.    1  and  simple  trigonometry,   we  write: 

SlN(cp-h&)=    -r-r 7^ (3-4) 

and 

cos((p  +  p)=  — p  -, <3-5) 

Substituting  equations  2.  5,    2.  10,    3.  1,    3.  2,    3.  3,    3.  4,    and  3.  5  into 
equation  2.  4,   we  have  the  radial  equation: 


a 


*V=      I       +e/?'~'l   ?""'■  1  (3.6, 


and  the  tangential  equation 
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W&     l/ZT^n?17        )  (3.7) 

We  get  some  simplification  in  equation  3.  6  by  substituting  for      & 
from  equation  3.7: 


V 


U      +  i,1  lzU*  r  (3.8) 

For  the  perturbation  solution,    we  take     $0     constant  and      U0      as 
given  by  equation  2.  12  to  write  the  assumed  solutions: 

U=U°    +    £U,-h£ZUz+  e7Uy  +-    •  •  •  (3#9) 

and 

■h  =  -U0  +  el,   +  ezK  +  e*l3  +  ■ '  •  {3.  io) 

where     U,  =  Uz -   •  '  •  -  #,  =  4\x=  •  -  •  =r  O       at     e,     ,    the  angle  where 
the  thrust  for  a  particular  cycle  is  started.     Substituting  equations 
3.  9  and  3.  10  into  equations  3.  7  and  3.  8,   we  obtain  the  first  perturba- 
tion differential  equations: 


j  /    __      -U^siAKp    +  U0   cos  <f> 
and 


(3.11) 


"i     ^U,--~T     +     «    '                 5-W^  (3.12) 

For  a  particular  cycle,   where  the  thrust  is  started  at  angle  0,      and 
G  -  6,   ■—   Z7T      ,    the  solutions  may  be  written: 

-&>)  =  J  <m^  (3.13) 
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and 

e  V  __, 

U/fi)^//i(e-J-|ji|((^         i*u*     StAj^l  Jf  (3ol4) 

For  a  unit  mass,   the  energy  is  given  by 

E    —      kinetic  energy  +  potential  energy 

Ro 

z. 


=  *.&['   -     ^r] 


2  ft*    J  <3-15> 

by  using  equations  2.  5,    2.  6,    2.  7,    2„  8,    and  2.  12,    and  where  the  po- 
tential energy  is  considered  zero  at     n   -  fc0    .     By  substituting  equa- 
tions 3.  9  and  3.  10  into  equation  3.  15,   we  obtain  the  first  perturba- 
tion energy  change: 

aJL  =  io(uzofu:i)i^(^u0-i)U  +.tfou:u;  0.16) 

The  perigee  distance  ratio  is  obtained  by  evaluating  equation 
2.  12  at        6  =  ©f 


R0  L^< 


'  o        


T~  (3.17) 


"We  could  require  that  this  quantity  remain  constant  during  the  thrust 
cycle,    as  done  by  Casey  in  reference  (2),   but  we  should  have  a  gain 
in  the  energy  imparted  if  we  require  that  the  perigee  distance  ratio 
be  equal  before  and  after  the  cycle.     We  shall  pursue  this  latter 
method  of  perigee  restriction  and  will  demonstrate  a  slight  gain  of 
energy  over  the  method  of  continuous  perigee  restriction. 
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The  eccentricity  may  be  evaluated  before  and  after  the  thrust 
cycle: 


e  =  \UuK'if  +  uk  <3-18> 


by  using  equation  2.  12.     Then,    substituting  equations  3.  9>    3.  10,    and 
3.  18  into  equation  3.  17  and  equating  the  results  before  and  after  the 
thrust  cycle,    we  obtain  the  first  perturbation  perigee  restriction: 

(3.19) 

evaluated  at      ©-     ,   where      O^       is  the  angle  at  the  end  of  the  thrust 
cycle. 
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IV.     OPTIMIZATION  OF  THE  RESTRICTED 
PERIGEE  CYCLE 

To  find  the  thrust  angle  which  will  impart  maximum  energy 
during  the  thrust  cycle  and  still  satisfy  the  perigee  restriction  given 
by  equation  3„  19,   we  shall  take  a  simple  variation,    using  the  La- 
grange multiplier,    equal  to  zero.     We  write  the  non-dimensional  en- 
ergy,   equation  3.  16,   plus  the  perigee  restriction  constraint,    equation 
3.  19,    with  the  Lagrange  multiplier,     A  ,    where      A    is  a  constant  to 
be  determined,    and,    assuming      6  =  constant  during  the  cycle,   take  a 
variation  on     <j>     equal  to  zero: 

For  convenient  notation,    we  rewrite  equations  3.  16  and  3.  19  as 

r~=  Afi,  +  Bu,  +cu;  (4.2) 

and 

PR  =  0  =   A%  +  Bu,  +  Cuf  {4a3} 

where     A    ,      A     ,    13   ,    and      C     are  shown  in  equations  3.  16  and  3.  19. 
Here,   we  note, by  talcing  the  difference  of  equations  4.  2  and  4.  3,    that 
AE 


£f° 


(=    (A-A')  t,(^)  (4.4) 

is  the  energy  change  for  the  thrust  cycle.     Substituting  equations  4.  2 
and  4.  3  into  4.  1,    we  have: 

(A +  !£')$&,  +U  +  D[BSUi  +  CSUi']  =  °  (4.5) 

Before  taking  the  variation,    we  must  first  integrate  equation  3.  14  by 
parts  to  obtain  the  more  useful  form: 
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u^)4^^-^-^^s/A/^^+^uCo5^-^w]^(y';j^  (4.6) 


Then,    taking  a  derivative  of  equation  4.  6  with  respect  to       Q 


& 


l.j.,,'1' 


& 


•J  n0    Mo 


|-,fs^^)<>^ 


(4.7) 


Using  equation  3.  11,   we  note  that: 


Ofc 


0, 


s</V=^V^^  = 


e, 


e, 


~  ul  Cos  (f>    -    Uo  SiA/^ 

i.uo3 1jl^[7uF~~~ 


$</>Jf  (4.8) 


Now,    taking  the  variation  of  equations  4.  6  and  4.  7,    using  4.  8,   we 
substitute  into  equation  4.  5  to  obtain: 

(4.9) 


For  equation  4.  9  to  be  zero,    we  must  require  that  the  integrand  be 
zero,    so  that 
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2  fZ- 


Substituting  the  expressions  for     /4     ,      J\'     ,    o      ,    and       C     from 
equations  3.  16,    3„  19,    4.  2,    and  4.  3,    using  equation  2.  12,    and  taking 


Op  —  O     »   equation  4.  10  becomes 

to.  n  <p  =  -  -ii*L£-  r — ri^ite')Lr — _ 


(4.11) 


where      A     is  a  constant  to  be  determined.     It  will  be  shown  that  the 
shift  of      ©p      is  negligible  for  the  symmetrical  thrust  cycle,    so  that 
a  choice  of  zero  for     Op       presents  no  problem  in  going  from  one 
thrust  cycle  to  the  next.     We  now  have  the  thrust  direction  program 
given  by  equation  4.  11. 

Substituting  the  expressions  for     A       and       A        and  equations 
2.  12,   3.  11,    and  4.  11  into  4.  4,   we  have: 
A  E        _ 

(4.12) 
Also,    substituting  equations  3.  11,    4.6,    4.7,    and  4.  11  into  3.  19,    and 
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using  2.  12,   we  have: 

PR-o  - 


aW 


(4.13) 


From  equation  4.  13  we  determine  the  value  of  the  constant       A<20 
for  any  particular  cycle,    since     n       always  occurs  with      S0       in  the 
two  equations  4.  12  and  4.  13. 

Since  it  is  not  possible  to  integrate  the  above  equations  exact- 
ly,  we  must  use  numerical  methods  to  carry  the  investigation  to 
large  angles.      This  would  be  fairly  simple  for  one  cycle,   but  when 
one  considers  a  large  number  of  cycles,    a  computer  program  appears 
to  be  mandatory.     We  shall  not  pursue  this  route,   but  instead,   we 
shall  present  an  approximate  evaluation  which  is  useful  to  angles  of 
one  radian. 
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V.  APPROXIMATE  EVALUATION  FOR  ONE  CYCLE 


To  make  an  approximate  evaluation  of  the  optimum  restricted 

perigee  equations  for  one  cycle,   we  shall  make  series  expansions  and 

integrate  term  by  term.     First,    using  equation  4.  13,   we  may  rewrite 

4.  1 2  as : 

{i+ef+te.coii'HH-e.cosi'fli+eZi-ze.cost'  +ZXe0(cos¥'- 1 )]    J^p 


*  6, 


DENOMINATOR 


(5.1) 


since  PR-  O  by  the  appropriate  choice  of        AS0       and  where 

DENOMINATOR    is  the  denominator  of  the  integrand  in  both  equations 
4.  12  and  4.  13. 

Now  consider  the  function 


G  =\  §*f 


(5.2) 


where 


h 


lj(i*e.coiV)t[LH-e*+ze.cosH')b  +Z7ieJcoiy-i)]2~  +  fef(n-e.?si*/zY 


( i  +e0toi  vo7  ]|  \+e?  +ie0cos<r 


(5.3) 


We  notice  that 


o2 


e7 


-L  l£ 


y-i 


9* 


JV 


b=i 


e, 


t=/ 


has  an  integrand  identical  with  equation  4.  13,    and 


9b 


2_ 

Db 


5H.51i.if. 


t  = 


t  =  l  ©. 


t*i 


(5.4) 


(5.5) 
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has  an  integrand  identical  with  the  non-vanishing  term  of  equation 
5.1,    since         PR-0      .      Thus,    if  we  perform  the  integration  of  equa- 
tion 5.2,    we  may  evaluate        /\eo        and  the  increase  in  energy  by 
taking  the  derivatives  indicated  in  equations  5.  4  and  5.  5. 

Expanding  in  series,   we  collect  the  coefficients  of  the  various 
terms  of     <£     and  integrate.     Noting  that  the  result  is  symmetrical 
about      Q  =  O        ,   we  write 

^KD     Cl^ep)L     z      Q+eJ-G     ^Q+Cflio     Ci+e.)4 10,090^  '"]      (5.6) 
where 

N,=  2e/te6-atr  +  (r2  (5.7) 


N2=  \Q>ef+  lde*-e?-e0  -22 cr<?f -32  <ree  +a<r+  ao  <r2ef  + 

^14^0  -4Taf  J2r3-3f4  (5.8) 


N 3  =  S++  ef  +  1036 ef+ 2 4Z e+-? t> 2 e* - 1  e ef  1- 2 e <,  -3+<?  re?-24iL  re?- 
-  H4C  re?-!-?**  re0  -  <f<r  +  /z  tz  rza?+  3£3SvzeI  -t28oz  <rze0z   - 
-3ZZTzeo  +•  32<rz  +  u+t>  (T3e*  +2&4o(r7e0-3oosr?  -^zor^e^  - 


(5.9) 


and  for  convenience  we  have  used 

(T  =     — 


b  (5.  10) 

Now,   taking  the  derivative  of  equation  5.  6  as  indicated  in  equation  5.  4 
equal  to  zero,   we  satisfy  the  perigee  restriction  of  equation  3.  19  and 
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solve  for     A  Q0    : 

Here,   we  note  that  equation  5.  1 1  would  be  much  more  complicated  if 
a  non- symmetrical  thrust  program  were  used.     However,    since  we 
can  see  from  the  form  of  equation  5.  6  that  this  will  have  no  energy- 
advantage,    a  symmetrical  thrust  program  is  chosen  to  simplify  the 
results  as  much  as  possible. 

Similarly,    taking  the  derivative  of  equation  5.  6  as  indicated  in 
equation  5.  5,    and  using  equation  5.  11,   we  have  the  energy  for  one 
half  of  the  cycle: 

j_    ae    _:     \]    (  (i-ze.)  &l        (7-  ueo  +  i&e?)     ef 

y       —      . _ —       ■+■ -^- 


.  r 

zizoe0  ±5i3oef+3<3Z4-e?-Z5O4-el-\50GeZ-\\34-eo-\\\i       &i 


+ 


( \  +ea)G  so^oo  ) 

(5.12) 

To  determine  the  advantage  of  this  method  of  perigee  restriction 

over  that  used  by  Casey,    the  energy  equation  given  in  reference  (2) 

was  expanded  in  series  and  integrated  to  the  same  number  of  terms 

as  we  have  in  equation  5.  12.      The  two  energies  are  identical  to  the 

last  term,   but  if  we  compare  the  last  term  of  the  constant  perigee 

energy  expression: 

zizoei  +5ieoef+\\zse*-zoooel-z55oe}-+ttoe0-z33f       el 

+ ^_        (5.13) 

0+eo)e  50*00 

with  the  last  term  of  equation  5.  12,   we  see  that  for  values  of      60 
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between  zero  and  one,    the  constant  perigee  restriction  yields  less 
energy. 

To  find  the  changes  of  the  elliptic  orbit  parameters      U  ,  -^  , 
and    Q.    ,    we  note  from  equation  4.  4  that 

V*,  =  (TV)F0r  (5-14) 

where  the  energy  is  given  by  equation  5.  12.  Substituting  the  expres- 
sions for  A  and  A  and  neglecting  the  higher  perturbations  since 
€     is  small,   we  have 

4&=£V%>=  T-^-T*   ^4  (5-15) 

0  +  e.)    $oRo 

Now,   from  the  perigee  restriction,    we  may  write  equation  3.  17  before 
and  after  the  thrust  cycle: 

I  +  Qo  I  +■  eo  +  A  e 


from  which  we  solve  for  the  change  in  eccentricity: 


(5.16) 


using  equation  5.  15.  To  see  if  there  is  a  shift  of  perigee,  we  assume 
that  the  equation  for  U  ,  at  the  end  of  the  thrust  cycle,  has  the  form 
of  2.  12  : 


U  = 


[I  +(ec-hAe)CDs(&z-*eP)] 


Z  (5.  18) 


where  we  have  provided  for  a  shift  of  the  perigee  with  the  term     A&p 
Then,   using  equations  2.  12,    3.  9,    and  5.  18,   we  can  write 
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*A  >2. 


/  r+  pi.       \ 

u,*il ^tt^  *  (5-19) 


and 


/        l 


£ 


U,=  c     ^ — ~T^ ?  (5-2°) 


[4,  +  A -A.) 


Substituting  these  two  equations  into  4.  3,    with  the  expressions  for 
A      *     B     t    and      £■     ,    we  obtain  the  result 

{cos(Gx-Aep)-cosat)cosez-t-(s'*i(6*-aep)-S"u61.)s"<'62,  =  O   +    0(€)  (5.21) 

Thus       A  9f>  ■      is  at  most  of  order      £     „     Using  this  result,   we  may- 
evaluate  the  terms  of  order     £      in  equation  5.  21,   using  5.  15  and 
5.  17,    and  we  find  that  they  are  zero.      Thus,    there  is  no  shift  of  the 
perigee  within  the  magnitude  of  the  first  perturbation  as  would  be  ex- 
pected with  a  symmetrical  thrust  program. 

Before  proceeding  with  the  discontinuous  thrust  program,   we 
shall  obtain  the  approximate  solution  for  the  spiral  orbit  with  thrust 
tangent  to  the  velocity  vector.     For  this  case,      <p     is  zero,    and  equa- 
tion 3.  8  become  s 

From  equation  3.  7  we  have 

Taking  twice  of  equation  5.  23  equal  to  the  derivative  of  5.  22  and  ex- 
panding,  we  obtain 
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where  we  have  assumed  starting  from  a  circular  orbit  with      (  U   /LA) 

and    (LA  / U)    much  less  than  one.     Integrating  from     O     to     &    we 

have: 

e 

6 

To  determine  approximately  the  small  quantities  not  yet  integrated, 
let  us  take  for  the  first  approximate  solution 

Ul  =    Uf  -  4<L&  (5.26) 

which  is  merely  the  first  result  from  equation  5,  25.      Then,    taking 
derivatives  of  equation  5.  26,   we  have  approximately: 

u  ul  {5- 27) 


u"  4£2 

U=--^  (5.28) 

UUM=   -    ^f'  (5.29) 

Substituting  equations  5.  27,    5.  28,    and  5.  29  into  5.  25,   we  obtain  the 
result: 

UZ  ^    Ul    +lZ&.Z\-[jr jjzj  (5.30) 

where       U^,     is  given  by  equation  5.  26.     We  may  now  use  the  second 
term  of  equation  5.  30  to  choose  values  of     G.     and       U^  ,    assuming 
U0      is  given,   to  insure  that  our  results  lie  within  a  certain  accept- 
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able  error.     We  should  also  note  that    (  &/ U^)      must  be  much  less 
than  one  if  equation  5.  30  is  to  have  meaning. 

To  determine  the  time,    we  use  part  of  equation  2.  13  and  sub- 
stitute equations  5.  22  and  5.  30.      Then,    using  equations  5.  26  and  5.  28, 
we  obtain 

KATE'S    **'*=  \  (*  +  %)*  M     eJo 


r-^-$-T*l-k'  *kti  +  --')*» 


4      nl  U/*       u:u*         u 


D 


(5.31) 


Here,  we  see  that  for  specific  values  of  Uq_  and  £  ,  we  will  have 
a  larger  error  in  equation  5.  30  than  in  5.  31,  since  IAA  is  less  than 
U0  .  and       U0     is  less  than  one. 

To  determine  the  energy  input,   we  substitute  equations  5.  22, 
5.  27,    5,  28,    and  5.  30  into  3.  15,    and  take  the  difference  between  the 
final  and  initial  values: 


A  E 


I 


3  H 


where  again  the  error  may  be  estimated  by  the  second  term.     Of 
these  three  equations,    5.  30,    5.  31,    and  5.  32,    equation  5.  30  will  have 
the  largest  error  as      UL&.      becomes  small,    and  it  should  be  used  for 
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choosing  values  of    6      and        U^     . 

We  have  reached  the  point  where  we  must  decide  how  to  pro- 
gram the  cycles  for  the  restricted  perigee.      To  do  this,    we  shall 
choose  a  ratio  of    (R0//l)  -0.9 5    to  be  the  maximum  value  for     U     . 
For  the  case  of  the  earth,   this  corresponds  to  an  altitude  of  208  miles 
and  is  adequate  to  insure  negligible  atmospheric  effects. 

If  we  form  the  quotient  using  equations  2.  13  and  5.  12, 

specific  power    =  7?  =    -~^    -A£-  (5o  33) 

then  expanding     A~T"     in  series  and  performing  the  division 


(5.  34) 


This  equation  is  plotted  in  fig.    2  for  various  values  of       Q.  G     ,   with 
corresponding  values  of     (E  / Qo  Ro)         for     P0  —  0.3  b         t   where 
1  o         is  defined  as 

-ho 
Also  plotted  in  fig.    2  is  the  specific  power  for  the  spiral  orbit  for  the 

same  values  of     (E  /  <^o  K  b)  .     We  can  see  that  when  the  eccen- 

tricity is  small,    there  is  not  much  difference  between  the  two  methods 
near  the  perigee,   but  the  discontinuous  thrust  program  becomes  inef- 
ficient if  we  proceed  too  far  from  the  perigee.      This  is  because  the 
thrust  must  be  directed  away  from  the  velocity  vector  to  satisfy  the 
perigee  restriction.      Thus,   for  a  given  energy  increase  per  cycle 
when  the  eccentricity  is  small,     &     must  be  larger  with  discontinuous 
thrust,    since  the  thrust  is  not  on  for  the  full      Z.TT     radians.      There- 
fore,    6  Y       is  larger,    and  the  payload  and  structure  to  initial  mass 
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ratio  is  reduced,    as  seen  in  equation  2.  25.     If  we  manipulate  equa- 
tions 2.  13  and  3.  15,    we  can  obtain  the  time  for      2  77"    radians  of  trav- 
el in  terms  of  the  total  energy: 

3 


AT     WiTL'      $°k 


^^-Al"^  (5.36, 


which  holds  for  both  elliptic  and  circular  orbits.     Now,    if  the  mission 
is  to  be  accomplished,    both  methods  will  have  about  the  same  number 
of  cycles,   where  in  the  case  of  the  continuous  thrust  method,    cycle 
means     £  77"     radians  of  travel.      Thus,  we  want  to  impart  about  the 
same  amount  of  energy  during  each  corresponding  cycle. 

We  can  also  see  in  fig.    2  that  for  larger  values  of  eccentricity, 
say  0.  75,    the  specific  power  is  much  greater  for  the  elliptic  orbit. 
In  this  case,   for  a  given  energy  increase  per  cycle,       /\Y     for  the 
spiral  orbit  will  be  greater. 
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VI.     APPLICATION  TO  MULTIPLE  CYCLES 

To  find  the  total  energy  input,    total  time,    and     Y   ,    we  must 
make  a  summation   of  the  results  for  each  cycle.     However,    since  we 
are  assuming  a  large  number  of  cycles,    with  small  changes  for  each 
cycle,    we  may  write  integrals  for  the  summation  operation  and  ob- 
tain our  results  by  integration.      The  thrusting  time,     At      ,    is  taken 
as  constant  for  all  of  the  cycles,    and  the  variable  of  integration  is 

Je     . 

We  must  first  express  the  cutoff  angle,      <3,     ,   in  terms  of  Q. 
because  it  must  be  varied  slightly  to  keep  the  thrusting  time  constant. 
Expanding  equation  2.  13  in  series,   we  have,   for  the  symmetrical 
orbit: 


Iff."  ?*ft^r[        Q  +  e)         i0° 


(6.1) 


Then,   if  we  let 

we  can  write  approximately 

^=ftfi^f<-ie&2  +  f(7fi7^  (6-3) 

where     jl     is  the  initial  value  of       Qz     .     Care  must  be  used  in 
choosing     -fe    so  that       Gz     is  less  than  one  at    £  -  I        ;  otherwise, 
the  series  expansions  will  not  be  good  approximations. 
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Substituting  equation  6.  3  into  5.  17,   using  5.  12,    and  replacing 
the  finite  "delta"  quantities  with  differentials,   we  obtain 

The  term     al        is  in  effect  one,    that  is,    one  cycle,    and     (<$&/€:)     is 
the  small  change  of  eccentricity  during  that  one  cycle.      To  perform 
the  integration,    we  shall  take     £     as  constant  to  obtain: 

cn-^**^*-  zlo0fT7r M       (6-5) 

where      £2/       is  the  initial  eccentricity  and      <£?r     is  the  final,      N     is 
the  number  of  thrust  cycles,    and  the  total  thrusting  time  is  given  by: 

6TT  =  £  NaT  (6.6) 

We  find  the  energy  input  very  simply  by  replacing  the  finite  "delta" 
quantities  in  equation  5.  17  with  differentials,  using  5.35,  and  inte- 
grating to  obtain 


Ef  -  Et  =  %£  Po  (ef  -  eL ) 


Then,   using  equations  2.  22,    3.  1,    and  6.  6, 


(6.7) 


ya=    i^2NAf.    ^^e'j^NAT  (6.8) 


To  find  the  total  travel  time  for      ZTf     radians,   we  write  equation 
5.  36  in  dimensionless  form,    replacing  the  finite  "delta"  quantities 
with  differentials;   then  substituting  3.  15  and  5.  3  5,   we  have 
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d  l  a -JFTTT?*  (6-9) 

where  again,       di       is  in  effect  one.     Substituting  equation  6.  4  for 
d  t       and  neglecting  "end  effects"  in  the  integration,    we  have 

6Ts4r(lt'W  ^z>w"  vr^1    I   (6-10) 

The  "end  effects"  we  mention  are  the  initial  cutoff  angle  and  the  final 
cutoff  angle  not  being  equal  so  that  we  have  a  small  amount  of  time 
neglected  for  each  cycle  due  to  the  shift  of  the  cutoff  angle;  the  final 
energy  is  added  before  starting  the  final  coasting  orbit;  and  the  initial 
energy  is  added  during  the  latter  part  of  the  initial  coasting  orbit. 
This  might  be  clearer  if  we  note  that  the  integration  gives  us 

where      A  ~Te  ■        is  the  travel  time  for     2TT      radians  at  eccentricity 


(6.12) 


Q :    y  whereas  the  true  mission  time  would  be  given  by: 

where         a"7^-  is  the  thrust  time  for  the  initial  cycle, 

&      —    S?  .  A  & .  ,    and       d  Tj        is  the  cutoff  angle  shift 

time.     If  the  number  of  cycles  is  not  large,   then  equation  6.  10  should 
be  corrected  for  "end  effects"  by  adding 

r     =  &t.  (i-h  5^— f)  -&~%.  -  A-re  (6.13) 

Com:  i  L  9t    i-k'  '  + 
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where    7£      is  the  initial  cutoff  angle  given  by  equation  6„  2,       &z  r      is 
obtained  by  evaluating  equation  6.  3  at      €.  -  £  f       .      The  cutoff  angle 
shift  time  correction  is  an  average  value  obtained  by  using  the  first 
term  of  equation  6„  1.      The  other  corrections  follow  simply  by  com- 
paring equations  6„  11  and  6.  12. 
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VII.     RESTRICTED  PERIGEE  CYCLE  WITH 
BATTERY  STORAGE 

We  must  still  consider  the  possibility  of  different  values  of 
specific  power  plant  mass,     K      ,    for  the  two  systems,    since  we  have 
the  capability  of  storing  energy  with  the  discontinuous  method     while 
coasting  if  we  provide  batteries.      The  energy  generated  by  the  power 
supply  while  coasting  could  be  stored  in  batteries  and  then  used  in 
addition  to  the  power  supply  while  thrusting. 

To  study  this  possibility,    we  shall  first  assume  that  a  battery 
pack  is  included  with  the  power  supply  and  that  other  conditions  of 
the  thrust  program  are  unchanged;  later  the  battery  will  be  consid- 
ered as  part  of  the  payload.     We  shall  further  assume  that  during  the 
early  part  of  the  program,   when    6      is  small,   that  the  coasting  time 
is  not  sufficient  to  fully  charge  the  batteries,   but  as     €.     increases, 
a  point  is  reached  where  the  batteries  do  become  fully  charged,    and 
are  fully  recharged  during  each  subsequent  coasting  period.     We 
have  equation  2.  18,   evaluated  at  the  "burnout  condition": 

Mf-  .  f  ^k  Jlf7""^  *L  AT  „  ,, 

Ml  M„     Z    I       K        P  (7-1) 

where  we  have  multiplied  and  divided  by     [K/ M^J  ,    and  must  now 
consider    P     as  a  function  of  time.     When  the  batteries  are  not  fully 
charged,   the  power  is  given  by 

r~    |<       AT  <7-2> 

neglecting  battery  losses  while  charging,    where     Ms        is  the  mass  of 
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Mt       '         M„     £    I       K        P  (7,1) 

where  we  have  multiplied  and  divided  by     (K/M^y  ,    and  must  now 
consider    P     as  a  function  of  time.     When  the  batteries  are  not  fully 
charged,   the  power  is  given  by 

P 


M<3         ATCWU- 


K         AT  (7-2> 

neglecting  battery  losses  while  charging,   where     M$       is  the  mass  of 
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the  electrical  generating  unit,    K     is  the  specific  mass  of  the  gener- 
ating unit,       A~f     is  the  constant  thrusting  time  per  cycle,    and 
AT        i       i-s  the  travel  time  for     2  7T      radians  given  by  equation 
5.  3  6.      The  batteries  will  be  fully  charged  when 


M  batt.  N, 


K  uB<+  K 


Korbit     -^)  t7.3) 


batt 


or 


>     ,_,      Mbatt.  K 


AT_v...^    ^T+  — (7.4) 

batt. 


M*  |< 


where    M,  and     K,     _       are  the  mass  and  the  specific  mass  of 

batt.  batt. 

the  battery,    respectively.     From  equation  7.  1,   we  define       Y      ,    cor- 
responding to  2.  22: 


/T.  ..  ,  T 


(7.5) 


But  since     6      is  constant,    and       H  uv    -     M<j   +-  Mba.th.  '' 

o 

where  -f  ^)         is  introduced  to  remind  us  that  the  integrand  is 

zero  when  not  thrusting.     Using  equation  7.  4  in  7.  2  ,    we  can  write 
the  maximum  power  available  when  the  batteries  are  fully  charged: 

*  max-    K    (}  M$       Ku+f.    AT   I  (7o7) 
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Using  equations  7.  2  and  7.  7  in  7.  6  : 
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,2    e  0 


r  J^  i  r. 
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J-i 


(7.8) 
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where        t      .         is  the  time  at  which  the  batteries  become  fully 
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charged.     Now,   the  integrals  in  equation  7.  8  may  be  treated  as  sum- 
mations since  the  thrust  is  on  for  the  constant  time      AT"     during 
each  cycle: 


K 

&     Ku-tr. 


(7.9) 


where       J\|      .,       is  the  cycle  number  after  which  the  batteries  are 
'^crit. 

fully  charged.     Using  equations  6.  9  and  7.  4,   we  can  solve  for  the  ec- 
centricity   when  the  batteries  become  fully  charged: 

Z1T 


ec~- 


f? 


AT  +    Mb^' 


(7.10) 


Then,   using  equations  6.  4  and  6.  9 
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Integrating  this  equation  and  determining       N  -  IV      it      from  equation 
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6.  5,    using  7.  10,    equation  7.  9  becomes: 

2    r 


y*  =  f  u^%^^\4^-^-m^-^ 
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where 
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I  +■  H   =    /    + 
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and 
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(7.  14) 


Now  if  we  form  the  ratio  of  equation  2.  22  for  the  spiral  orbit  and 
equation  7.  12,    multiplied  by  the  factor         J   f-  ^/         ,   we  have 
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FRYz=(n-/*) 
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When      FRy^         is  less  than  one,    the  spiral  method  of  energy  addi- 
tion must  have  a     Y     less  than  the     Y     for  the  restricted  perigee 
method,    and  by  equation  2.  25  the  spiral  method  must  have  a  greater 
payload.     If       FRyZ         is  greater  than  one,    then  there  may  be  a  pos- 
sibility of  carrying  a  greater  payload  by  using  the  restricted  perigee 
method.      This  will  depend  on  the  ratio  of        K/Kk   ++•       that    can  be 
obtained. 
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VIII.     COMPARISON  OF  RESULTS 

Now  to  compare  the  results  we  have  obtained,    equation  5„  31, 
the  spiral  thrust  time,    is  plotted  against  the  energy  added,    equation 
5.  3  2,    in  fig.    3.     Using  a  value  of         K  -  0.  8    ,   which  corresponds  to 
a  maximum  6^  =  0.  944    by  equation  6.  3,    the  restricted  perigee 

thrust  time    is  computed  using  equations  6.  2,    6.  5,    and  6.  6  for  values 
of  eccentricity  from  zero  to  0.  95,   with  F^   =  0.  95  .      The  corre- 

sponding increase  of  energy  is  computed  using  equation  6.  7,    and 
these  results  are  plotted  in  fig.    3.     Since  the  total  time  must  be  equal 
for  either  method  at  a  specified  energy  addition,   we  may  compute 
equations  5.  31  and  6.  10  at  specific  energy  levels  and  take  their  ratio 
to  plot     (€..  /£s)  in  fig.    3.      This  is  the  ratio  of  acceleration  for 

the  discontinuous  to  the  spiral  that  is  necessary  to  accomplish  the 
mission. 

It  seems  clear  from  fig.    3,   in  view  of  the  acceleration  ratios 
required  for  the  mission  and  equations  2.  13,    2.  22,    and  6.  8,    that    T 
will  be  greater  for  the  discontinuous  method  if     K    is  the  same  for 
both  methods  and  starting  from  the  circular  orbit  with  F^   -    0.  95  . 

For  these  conditions,    the  discontinuous  thrust  method  will  provide 
less  payload.     Further,   if  we  decrease     7t     ,    the  initial  cutoff  angle, 
we  see  from  equation  6.  10  that     £    must  be  increased  for  the  same 
total  time;  thus,        T        will  be  even  larger  for  the  discontinuous 
method,    resulting  in  even  less  payload. 

Using  equation  7.  10  and  fig.    3,   we  have  evaluated  equation 
7.  15  for  three  different  final  energies,    expressed  as  eccentricity,   for 
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various  values  of    o<  .      These  results  are  plotted  in  fig.    4. 

From  reference  (3),   we  find  that  a  typical  value  of     K     would 

be  about  150  pounds  per  kilowatt,,      This  value  will  decrease  as  space 

power  plant  technology  progresses.     In  reference  (3),    we  also  find 

Kl    tt       of  7  5  pounds  per  kilowatt-hour  for  nickel-cadmium  batteries, 

which  have  very  reliable  charge  and  discharge  characteristics,    and 

14  pounds  per  kilowatt-hour  for  silver-zinc  batteries,   which  are  not 

as  reliable  for  a  large  number  of  recharging  cycles.     From  equation 

7.  3,    we  see  that  both  values  of       K-,     -^       must  be  divided  by 

batt.  ' 

3600(<3o/Ro)  to  have  the  proper  non-dimensional  time  units. 

Thus,  we  have  the  values  for       W,     .,       of  16.  8  for  nickel-cadmium 

,vbatt0 

and  3.  13  for  silver-zinc. 

We  cannot  improve  over  the  spiral  thrust  program  with  either 

battery  supply  if  the  final  energy  must  correspond  to  an  eccentricity 

of  0.  95,    since  in  fig.    4  the  maximum  (     RVZ     is   !•  00&  • 

0 

Using  equations  5.  35  and  6.  9,    we  compute  the  non-dimension- 
al thrust  time  for  one  cycle,    corresponding  to  an  initial  cutoff  angle 
of  0.  8  and        eo  =  O         : 

r      Cmo.*)£r=     -^-j=    i.728  (8.1) 

(  ZT      J1         (035)  /z- 

Using  this  and  equation  7.  14,   we  compute  the  values  of     yU    for  the 
nickel -cadmium  battery,   for    <=<  =  7.  5    and       o<  =  10.  0,    respectively: 

y  =   l.+f     }        p^  1.33  (8.2) 

and  for  the  silver-zinc  battery: 

fA  =  0.27     ,        y  -0.3£  (8.3) 
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From  fig.    4,   we  read  FRy2  =  l.  21    for  the  final  eccentricity  of 

0.  4  and       ex  =s  7„  5  0     But  from  equations  7.  15,    8.  2,    and  8.  3,    the 
values  of    ^      for  either  battery  will  require         oi  to  be  greater 

than         T$  ;  again,   we  find  no  advantage  in  the  restricted  perigee 

method  of  energy  addition.     However,   for  the  final  eccentricity  of  0.  7 
we  read  from  fig.    4  FRyZ     of  1.  29  and  1.  39  for    o<'s  of  7.  5  and 

10.0,    respectively.     From  equations  7.  15  and  8.3,   we  compute  for 
this  condition: 

for  the  two  values  of     c<  ,    respectively.     From  equation  2.  25,   we  can 
see  that  in  this  case  the  restricted  perigee  method  has  a  slightly 
larger  payload  than  the  spiral  method  for  tins  particular  final  energy. 
To  see  this  result  more  clearly,   we  form  the  ratio  of  equation  2.  25 
for  the  two  methods  using  8.  4: 


2 


(8.5) 


for      ex!  —  10°  0  •     Now,    clearly,    the  numerator  of  8.  5  is  greater  than 
the  denominator,    so  the  ratio  must  be  greater  than  one. 

Although  we  have  found  that  the  restricted  perigee  method  of 
energy  addition,    using  silver-zinc  batteries,   yields  a  slightly  larger 
payload  than  the  spiral  method  for  certain  values  of  final  energy,   we 
still  have  the  problem  of  battery  reliability.      There  is  no  question  of 
the  recharging  reliability  of  the  nickel -cadmium  battery,   but  we 
found  the  specific  mass,      Ki,    f,     >   to  T°e  to°  bigh  to  be  useful.      The 
silver-zinc  battery  could  be  used  for  a  slight  payload  advantage,   be- 
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cause  of  its  lower  specific  mass,    but  its  recharging  reliability  must 
first  be  improved. 

We  must  finally  compare  the  iwo  methods  of  energy  additio  1 
when  reliable,    rechargeable  batteries  are  part  of  the  payload  but  are 
available  to  be  used  in  the  thrust  program.,     In  this  case,    LA  =  O      a-nd 
we  see  immediately  from  fig.    4  that  there  are  certain  values  of 
and  final  energy  where  the  restricted  perigee  has  a  decided  advantage 
over  the  spiral  method. 

As  an  example,    let  us  determine  the  advantage  where  a 
nickel-cadmium  battery,       K-.     ++     =   16.  8  ,    is  part  of  the  payload. 
Using  equations    7.  14  and  8.  1,    we  compute  values  for    (^uff.  / ^%) 
for     ex  =  7.  5    and      <x  =  10.  0,    respectively: 

=     1.45  _Mba£_     =    !.93  (8.6) 


Even  though  the  ratio  of  battery  to  power  supply  mass  is  rather 
large,    it  is  not  unrealistic  since  we  could  easily  have  a  satellite 
where  both  battery  and  power  supply  are  but  a  small  fraction  of  the 
initial  mass.      From  fig.    4  we  read  PR  v2  =  1«  21  for  the  final 

eccentricity  of  0.  4  and     c<   =  7.  5  ,    and  since        u  -  O        '• 

Forming  the  ratio  of  equation  2.  25  for  the  two  methods,   using  8.  7  : 

(J±J^lk=     (I  -  ^  )Z  (3.8) 

(M,+Ms)s        (1-l.lOYj)2 

Similarly,    from  fig.    4  we  read         FRyz       of  1.  29  and  1.  39  for  the 
final  eccentricity  of    0.7  and      ex's  of  7.  5  and  10.  0,    respectively. 
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Since      fJ  =  0      ,    we  have  immediately: 


-£-=l.l3C       ,     -f 


=  I.I3G       ,     -£■  =  1.179  (8.9! 


for  the  two  values  of     ex.    respectively.     Again  forming  the  ratio  of 
equation  2.  25  for  the  two  methods,    using  8.  9>    we  have 

(M,  +  Ms)s         (|-i.|36Xi)Z     '(M,  <-M,)s-  0-l.'79%l)8  {8- 10) 

for  the  two  values  of      oi  =  7.  5    and      o<  =  10.  0,    respectively.      The 
ratios  in  equations  8.  8  and  8.  10  will  depend  on  the  value  determined 
for       V.      ,    computed  from  equation  7.  12,    but  they  do  show  a  definite 
advantage  for  the  restricted  perigee  method  of  energy  addition. 

For  example,    suppose  the  value  for      Yj       is  0.  25,   that  is, 
let  25  per  cent  of      MQ     be  propellant.     Substituting  in  equation  8.  10, 
we  compute  ratios  of  1.  097  and  1.  131  for  the  two  values  of      o<    ,    re- 
spectively.    Now,    this  is  about  a  10  per  cent  gain  in  payload  and 
structure  mass  by  using  the  discontinuous  rather  than  the  continuous 
method  of  thrust  programming,    where  the  batteries  are  part  of  the 
payload. 

It  seems  rather  surprising  that  we  have  found  no  advantage 
for  the  restricted  perigee  method  of  energy  addition  for  the  final  ec- 
centricity of  0.  95,    where  we  found  a  definite  advantage  for  the  final 
eccentricity  of  0.  7.      There  are  two  reasons  for  this  result.      First, 
we  can  see  from  equation  6.  9  that  the  travel  time  for     2.  IT    radians 
increases  as  the  eccentricity  increases.      Thus,    if  we  thrust  for  the 
same  amount  of  time  during  each  cycle,    there  is  a  greater  percentage 
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of  coasting  time  as  the  eccentricity  increases.      This  requires  an  in- 
crease of        (^■d/^-s)  '    as  shown  in  fig.    3,    to  accomplish  the 
mission,    and  in  turn,    we  find  a  decrease  in  the  ratio  of      (Ys  /  Y^  ) 
so  that  there  is  no  advantage  when      B^       reaches  0.  95  .      Second,    by- 
restricting  the  cutoff  angle  to  one  radian  or  less,    we  do  not  take  full 
advantage  of  the  greater  velocities  provided  by  the  elliptic  orbit  as 
the  eccentricity  approaches  one.      This  can  be  seen  clearly  in  fig.    2 
where, at  a  cutoff  angle  of  one  radian  and      €?    =  0.  95  ,    the  specific 
power  of  the  restricted  perigee  is  many  times  greater  than  the  spiral 
method  of  energy  addition.      To  obtain  the  full  advantage  of  this  much 
greater  specific  power,    we  must  go  to  larger  cutoff  angles. 

A  study  of  the  restricted  perigee  method  of  energy  addition 
at  larger  cutoff  angles  would  require  the  use  of  numerical  methods. 
This  would  be  rather  involved,    using  equation  4.  13  to  evaluate    A  <S 
for  various  values  of        Q?      anc^      ^     >    then  using  this  result  for  a 
corresponding  evaluation  of  the  change  in  energy,    equation  4.  12. 
After  tabulating  these  "cycle"  results,    summations  of  time  and  energy 
could  be  made  with  a  corresponding  calculation  of       Y        for  various 
values  of    °<  .     By  making  such  a  numerical  study,    we  can  clearly 
define  the  area  where  the  restricted  perigee  method,    with  batteries, 
has  an  advantage  over  the  spiral  method  of  energy  addition.     We  will 
also  be  able  to  determine  if  there  is  any  area  where  an  advantage  ex- 
ists without  the  use  of  batteries. 
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IX.     CONCLUSION 

We  have  found  that  the  discontinuous  method  of  energy  addi- 
tion,   using  a  restricted  perigee,    has  no  payload  and  structure  mass 
ratio  advantage  over  the  spiral  method  when  a  specified  energy  is 
imparted  in  a  specified  time  and  the  cutoff  angle  is  restricted  to  one 
radian  or  less.      There  is  the  possibility  of  a  very  small  advantage, 
for  a  certain  range  of  final  energies,   if  silver-zinc  batteries  are  in- 
cluded in  the  power  supply  to  store  energy  during  the  coasting  periods, 
providing  their  reliability  can  be  improved  to  an  acceptable  level.     On 
the  other  hand,    a  moderate  decrease  of  the  power-plant  specific  mass 
will  rule  out  any  advantage  of  including  batteries  in  the  power  supply. 

In  the  case  where  batteries  are  included  in  the  payload,    say, 
to  periodically  energize  a  transmitter  after  the  satellite  has  been 
placed  in  final  orbit,   we  can  definitely  increase  the  payload  and 
structure  mass  ratio  for  certain  values  of  final  energy  by  using  these 
batteries  in  the  restricted  perigee  thrust  program.      This  does  not  ap- 
ply in  general,    obviously,    since  the  mass  fraction  of  the  batteries  re- 
quired in  the  payload  must  be  compatible  with  the  mass  fraction  of 
batteries  required  for  the  thrust  program  to  accomplish  the  mission. 
When  we  do  have  a  compatible  battery  requirement,    then  we  can  in- 
crease the  payload  and  structure  mass  ratio  by  using  the  restricted 
perigee  method  of  energy  addition. 

Since  we  are  restricting  our  investigation  to  low  thrust  de- 
vices,   a  numerical  study  of  the  restricted  perigee  method  of  energy 
addition  at  cutoff  angles  greater  than  one  radian  could  be  made  usino- 
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the  perturbation  equations  developed  earlier.     Such  a  study  would 
follow  the  procedure  briefly  outlined  in  the  last  paragraph  of  the  pre- 
ceding section,    and  is  necessary  to  clearly  define  the  area  where  the 
restricted  perigee  method  can  be  used  to  advantage. 
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